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We compute the two-point angular correlation function w(9) for a sample of ~ 1700 
galaxies to a magnitude-limit equivalent to R ~ 29.5 using a catalog derived from the 
Hubble Deep Field images. A non zero value of w(9) is measured down to R = 29.0. The 
amplitude of w(9) at the bright magnitude limit (R ~ 26) is consistent with previous 
ground-based observations. At fainter magnitudes the clustering amplitude continues to 
decrease but at a slower rate than that predicted by the power law w(l") oc io-°- 27R 
observed for shallower samples. The observed w(9) over the magnitude range 20 < 
R < 29 is consistent with linear evolution of the clustering of a galaxy population 
which at present has a correlation length r of about 4 h~ l Mpc, close to that of local 
IRAS galaxies. We also investigate the impact that magnification bias induced by weak 
gravitational lensing may have on our results. Although the observed amplitude of 
w(9) can differ from the true amplitude by up to 30%, this effect is not large enough 
to affect our conclusions. Finally, by using a color-selected sample, we examine whether 
the expected effects of magnification bias can be used for an independent determination 
of cosmological parameters in deep images. We conclude that the amplitude of the effect 
can be large and in some cases even produce an upturn of the amplitude of the correlation 
with limiting magnitude. However, we find that it is not possible to detect the effects of 
magnification bias on w(9) from images alone. If redshift information becomes available, 
it is possible to measure the effects of magnification bias directly and thus constrain the 
density parameter fl and the bias factor b. 



Subject headings: cosmology: observations, gravitational lensing, large-scale structure 
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1. Introduction 

An important constraint on the formation and evo- 
lution of structures in the Universe is the three di- 
mensional two-point correlation function £(r) 
function of redshift z. Unfortunately, redshift sur- 
veys to measure the intrinsic clustering properties of 
faint galaxies (R <^25.5), presumably at high red- 
shifts (z ^ 1), are difficult even with the new genera- 
tion of large aperture telescopes. Such direct mea- 
surements of £(r) for z ~ 1 are only now becom- 
ing available (eg. Cole et al. 1994, Le Fevre et al. 
1996). Therefore, to study the clustering properties 
of faint galaxies, one must for now rely on studies 
of the angular two-point correlation function w(9) . 
Some constraints on the redshift dependence of £(r) 
can be obtained by investigating the dependence of 
the amplitude of w{9) on the limiting magnitude. 

Recent studies of w(9) have pushed the limiting 
magnitude to ever fainter flux levels. Current ob- 
servational limits reach R = 26 (eg. Brainerd, Smail 
& Mould 1995, hereafter BSM). There is a general 
agreement that the amplitude of w(9) decreases more 
rapidly with limiting magnitude than expected from 
a redshift distribution N(z) as predicted by "no- 
evolution models" and linear evolution of the clus- 
tering. However, the interpretation of these observa- 
tions depends on the assumed model for the redshift 
distribution and clustering evolution, which are both 
poorly constrained by current data (eg. Glazebrook 
et al. 1995, Le Fevre et al. 1996). Some authors have 
argued that good agreement with the data can be ob- 
tained with models which assume modest clustering 
evolution of locally observed low surface brightness 
galaxies (BSM). 

Extending the analysis to fainter magnitudes is of 
great interest in order to impose more stringent con- 
straints on the epoch of formation of structures. How- 
ever, as first pointed out by Villumsen (1996, here- 
after V96), the interpretation of such data needs to 
take into account the effect of magnification bias in- 
duced by weak gravitational lensing, which may affect 
the measurement of w(9) at faint magnitudes. As dis- 
cussed by V96, this effect is expected to be important 
for samples with a median redshift 1 and it may 
therefore affect the analysis of very deep galaxy sam- 
ples such as those extracted from the Hubble Deep 
Field (hereafter HDF, Williams et al. 1996). The ef- 
fect should be most evident in samples which pref- 
erentially include "red" galaxies and therefore have 



shallow number counts slopes (Broadhurst 1996). A 
possible signature of the effect would be an upturn of 
the correlation amplitude with the median redshift of 
the sample which should correlate with the magnitude 
limit. 

The detection of magnification bias could be an im- 
portant tool to further constrain cosmological models. 
The amplitude of this effect is a measurement of the 
clustering of the mass and depends on the product 
f2o°8j where fio is the cosmological density param- 
eter and o\ is the variance of the moss fluctuations 
within a sphere 8 h^ 1 Mpc in radius. Therefore, the 
behavior of the clustering amplitude as a function of 
the limiting magnitude could provide a test on the 
value of the product flo&s- 

In this paper we use the HDF to investigate the be- 
havior of the amplitude of w(9) with magnitude reach- 
ing at least 3 magnitudes fainter than published data 
from ground-based observations. The faint magni- 
tude limit and the color information make these data 
ideal, except for the small angular coverage, to in- 
vestigate the behavior of w(9) at faint flux levels and 
the contribution of magnification bias to the observed 
clustering of faint galaxies. Previous work on galaxy 
clustering in the HDF field has been carried out by 
Colley et al. (1996) focusing on very small angular 
scales and discussing the possible existence of sub- 
galactic clumps at high redshift. Here, instead, we 
use the HDF data to investigate the evolution of £(r) 
as a function of redshift. 

In section 2 we describe the catalog used in the 
analysis. In section 3 we predict the correlation func- 
tion w(9) in the absence of magnification bias for sim- 
ple models of the redshift distribution and clustering 
evolution. In section 4, we compute the two-point 
angular correlation function for different magnitude- 
limited samples and determine the variation of the 
clustering amplitude as a function of the magnitude 
limit. Section 5 describes the theoretical calculation 
of the effects of magnification bias on w(9) , and com- 
pares the modified curves with the data. In section 6, 
the same analysis is done for a color-selected sample 
which should be more sensitive to the effects of the 
magnification bias. Our conclusions are summarized 
in section 7. 

2. The Galaxy Catalog 

Williams et al. (1996) presented a catalog of galax- 
ies extracted from the HDF images with FOCAS. An 
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alternative galaxy catalog has been used by Clements 
& Couch (1996). This catalog which was generated 
using the SExtractor program (Bertin & Arnouts, 
1995), was kindly provided to us by Couch (1996). 
As pointed out by Williams et al. , the FOCAS cat- 
alog finds in a significant number of cases several ob- 
jects where visual inspection of the images indicates 
that there is only a single galaxy. A comparison of 
the two catalogs with the HDF images leads us to be- 
lieve that Couch's catalog contains a smaller number 
of such cases. For this reason, we used that catalog 
for this work. A minimum object extraction area of 
30 pixels and a detection threshold of 1.3 u above the 
background were used. The magnitudes were com- 
puted with the zero points given by Holtzman et al. 
(1996). For this work, we used the catalogs extracted 
from the images taken with the F606W filter, which 
is similar to an R passband, and the F814W filter, 
which corresponds to the I band. Since only a very 
small number of stars is expected in the HDF field, 
we did not attempt to use the galaxy/star separa- 
tion parameter given by the SExtractor program and 
treated all detected objects as galaxies. In order to 
avoid edge problems, only galaxies with pixel coor- 
dinates in range 250 < x, y < 2050 were used. The 
total of 1732 galaxies detected in the F606W filter 
were used, out of which 1256 were detected in both 
the F606W and F814W filters. Hereafter, we refer to 
the two bands as R and I, respectively. 



where zq is approximately the median redshift (eg. 
Efstathiou et al. 1991) and Y is the Gamma function. 
Values for z were provided by Chariot (1996). The 
adopted redshift distribution predicts that for a mag- 
nitude limit of R = 28, 84% of the galaxies are at 
z > 1 and 29% of the galaxies at z > 2. Given the 
uncertainties, these estimates are in reasonable agree- 
ment with the estimated redshift distribution based 
on photometric redshifts from HDF (Mobasher et al. 
1996). It is important to note that the redshift distri- 
bution and hence the median redshifts are quite un- 
certain for faint, magnitude-limited samples such as 
the one considered here. This is probably the largest 
uncertainty in inferring the amplitude and evolution 
of £(r,z) from the present sample. 

With this parameterization, we can calculate 
w{9) using Limber's equation (eg. Efstathiou et al. 
1991) 
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Here H(z) is the Hubble constant as a function of red- 
shift, normalized so that H(z = 0) = 1, and x(z) = 

2 (l - (1 + z)~i) is the comoving angular diameter 
distance. Since the results for the observed correla- 
tion function depend only weakly on fio (BSM), we 
assume hereafter that Qq = 1. 



3. Predicted Correlation Function 

The method to derive the angular correlation func- 
tion from the intrinsic correlation function (neglecting 
magnification bias) for a given redshift distribution is 
well known, (eg. Peebles 1980). Following Efstathiou 
et al. (1991) we assume that the evolution of the in- 
trinsic correlation function is given by 
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when expressed in proper coordinates. A power index 
e = 0.8 corresponds to linear evolution of the correla- 
tion function, while e = corresponds to a correlation 
function constant in proper coordinates. Here ro is 
the present day correlation length. 
We adopt the redshift distribution 
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4. Correlation Function from the HDF Cata- 
log 

We have extracted from our catalogue eight R 
magnitude-limited samples with magnitude limits 
ranging from R = 26.0 to 29.5 in 0.5 magnitude steps, 
discarding all galaxies brighter than R = 23. Al- 
though the samples as defined are not totally inde- 
pendent, they are nearly so because the sample size 
increases rapidly with limiting magnitude. 

The angular correlation function w(9) is estimated 
using the estimator described by Landy and Szalay 
(1993) and BSM 

, m DD -2DR + RR 
W{ ) = RR • ( } 

Here DD, DR, RR are the number of data-data, 
data-random, and random-random pairs at a given 
angular separation. 

For each magnitude limit we generated a random 
sample within the same region of the galaxy catalog 
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but with five times as many objects. The number of 
DR and RR pairs are scaled to the number of DD 
pairs. The angular correlation function w(0) for each 
magnitude limit was estimated from pairs of galaxies 
with angular separations in the range 2'' < < 80". 
The total number of pairs at a given separation was 
calculated by summing the number of pairs at the 
corresponding separation within each individual chip. 
Pairs across chip boundaries were not included in or- 
der to minimize additional uncertainties associated 
with chip-to-chip variations of the photometric zero- 
point. The finite number of galaxies in the random 
samples add an uncertainty to the observed correla- 
tion function. However, this uncertainty is far smaller 
than the uncertainty due to the finite number of real 
galaxies. A recalculation of w{0) with ten times as 
many random galaxies as real galaxies lead to the 
same result. As a consistency check, we have also es- 
timated w{9) from counts in cells and the results are 
similar to those shown here. 

Errors were estimated from both Poisson statistics 
and from 30 bootstrap resamples of the data (see eg. 
Barrow, Sonoda & Bhavsar 1984). Both error esti- 
mates agree well at large separations. Fits to the 
correlation function are not sensitive to the differ- 
ences between the error estimates at small scales. The 
errors estimated from the bootstrap resampling are 
consistent with the field-to-field variations of w{0) as 
measured for the three chips. 

Due to the small angular size of the chips, it is nec- 
essary to take into account the "integral constraint" . 
The background density of galaxies is estimated from 
the sample itself, forcing the integral of the correla- 
tion function over the survey area to be zero. Since 
the angular size of the survey area is small, w(0) is 
reduced by the amount 

C=^J J dShdShwifi), (5) 

the so-called integral constraint (BSM). Here f2 is the 
solid angle of the survey area on the sky. If we assume 
that w(0) is a power law 

w(0) = A 9~ 1+1 ■ 7 = 1.8, (6) 

then C = 0.071A for our survey geometry for mea- 
sured in arcseconds. 

Figure 1 shows the observed w(8) for the eight 
magnitude limits. We determined the amplitude of 
w{0) by fitting 

w{0) =A0-~< +1 -C- 7= 1.8, (7) 



which takes into account the integral constraint. 
These fits are shown as solid lines in figure 1. The 
data points at different angular separations are cor- 
related but this is ignored in the fit. The error bars 
represent 1-er Poisson errors. There are two possible 
systematic effects that may affect our w(0) estimates. 
The first is due to merging of images, especially at 
fainter magnitudes. This effect should not be impor- 
tant in our calculations as they will preferentially af- 
fect w(0) on very small angular scales which have lit- 
tle weight on our fits to the data. A second potential 
effect is due to the incompleteness of the galaxy sam- 
ple close to the magnitude limit of the HDF data. 
Although there is no detailed study of the incom- 
pleteness of the HDF galaxy sample, the behavior of 
the number counts down to R ~ 29 indicates that 
incompleteness should not be important at brighter 
magnitudes. Moreover, the conclusions presented be- 
low would not change if we considered only galaxies 
brighter than R ~ 28 which is 1 magnitude brighter 
than the magnitude at which we may reasonably ex- 
pect the onset of incompleteness effects. 

We have applied Kendall's r test to estimate the 
significance level at which the null hypothesis of zero 
signal can be rejected. We find that the significance 
of the detection is > 2er for all subsamples, except 
the brightest, presumably because the small number 
of galaxies, and the faintest, possibly because of the 
incompleteness of the sample. In Table 1 we summa- 
rize the observational results. The first five columns 
give the magnitude limits, the estimate of the median 
redshift, the number of galaxies in the sample, the 
amplitude of w(8) at 1", and the associated la error. 

Figure 2 shows the amplitudes of the correlation 
function as derived from the fits of figure 1 as a func- 
tion of magnitude limit, scaled to a separation of 
1 arcsec. The filled circles are the HDF data and the 
error bars are 1-a uncertainties from fitting Eq. (|7]) to 
the measured correlation function. The open circles 
are the data taken from BSM but converted to the 
amplitude at 1 arcsec. The figure includes data from 
Couch et ai. (1993), Efstathiou et ai. (1991), Roche 
et al (1993), and Stevenson et al (1985), in addition 
to the BSM data points. 

The results of the HDF data are consistent with 
those of BSM at bright magnitudes, despite the ad- 
mittedly large error bars. However, at fainter lim- 
its the clustering amplitude falls off slower than the 
power-law A(R) oc lO" 027 R proposed by BSM 
to fit the data at brighter magnitudes in the range 
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18 < R < 25. Although all HDF points lie well 
above this line there is no indication of an upturn at 
very faint magnitudes as previously claimed by dif- 
ferent authors in the B-band (eg. Neuschaefer et al. 
1991, Landy, Szalay, & Koo 1996). 

The curves in figure 2 are theoretical predictions 
of the amplitude of w(6) assuming a redshift distri- 
bution given by equation ^, using values for the me- 
dian redshift as listed in Table I. We have considered 
models with e = and 0.8, and ro = 2 to 5 h^ 1 
Mpc in 1 hr 1 Mpc steps. The dotted curves are the 
predictions for a local galaxy population with differ- 
ent correlation lengths r and e = 0, while the solid 
curves are the predictions for e = 0.8 . The theoretical 
curves predict that the correlation amplitude will de- 
crease slower than the power law observed at brighter 
magnitudes, consistent with our data. 

A model with linear evolution of the clustering am- 
plitude (e = 0.8) with a present day correlation length 
ro ~ 4 h" 1 Mpc fits the data remarkably well over a 
range of nine magnitudes, 20Jsi?Js29. This value of 
ro is very similar to that measured for IRAS galaxies 
(Fisher et al. 1994), consistent with the picture that 
faint HDF and present day IRAS galaxies are drawn 
from the same population of field galaxies. If this is 
true, it would suggest that IRAS galaxies were formed 
at redshifts z ^ 1.5. As pointed out above, these 
conclusions depend critically on the assumed redshift 
distribution and must await supporting evidence. 

This conclusion contrasts to that of BSM who ar- 
gue in favor of a very weakly clustered population 
with ro ~ 2 h~ x Mpc, presumably low-surface bright- 
ness galaxies, or the interpretation Landy, Szalay, & 
Koo (1996) for a low redshift population of very faint, 
weakly clustered blue galaxies. The reason for this 
discrepancy is the higher assumed median redshift, 
which is supported by the results of Mobasher et al. 
(1996). At faint magnitudes, the HDF data can also 
be fitted with a model with slow evolution i.e. e = 0, 
if ro = 2 to 3 h~ x Mpc. However, such a model would 
be inconsistent with the data for R^, 25, being too 
shallow at brighter magnitudes. 

Since we have assumed that the median redshift 
of galaxies in HDF is larger than unity, the effects of 
weak gravitational lensing can influence the observed 
correlation function of galaxies. This so called mag- 
nification bias could in principle affect our conclusion 
(V96). The importance of this effect for the HDF 
sample is evaluated in the next section. 



5. Effects of Magnification Bias 

V96 has shown that the magnification bias can 
have a significant influence on w(ff) provided that: 
1) the sample depth in redshift is zjSi 1; 2) the slope 
of the number counts s is significantly different from 
0.4; and 3) the quantity f2oO"8 is not much less than 
unity. We emphasize that ag is the rms mass density 
fluctuations on a scale of 8 /i" 1 Mpc. 

The magnification bias, as shown by Broadhurst, 
Taylor & Peacock (1995), will affect the number den- 
sity of galaxies at a given position on the sky. If the 
slope of the number counts is s and the magnification 
due to gravitational lensing from a matter overden- 
sity is fj,, then the observed number density N t, s of 
background objects will be different from the "true" 
number density N trU e of objects 

N obs = N true fi 2 - 5 "- 1 . (8) 

Note that fi = 1 corresponds to no magnification. As- 
suming there are no intrinsic correlation of galaxies, 
then apart from random fluctuations due to the finite 
number of galaxies, the number density of galaxies is 
constant across the sky. In the presence of mass den- 
sity fluctuations, gravitational lensing will change the 
number density of galaxies according to Eq. @ . The 
amplification will be a function of position in the sky, 
and therefore the number density of galaxies will be a 
function of position. This means that a non-zero cor- 
relation amplitude is observed even in the absence of 
intrinsic galaxy clustering. In the limit of weak clus- 
tering, i.e. \fj, — 1| <C 1, the magnification /i relates to 
the "convergence" k as /i = 1 + 2k (eg. Blandford et 
al. 1991, Broadhurst, Taylor & Peacock 1995). The 
quantity k is a dimensionless measure of the surface 
mass overdensity. If the intrinsic relative overdensity 
is A/i we find that 

N obs = A (1 + A/i + (5s - 2)«) . (9) 
Here Nq is the average number density of galaxies so 
that Ntrue — N Q (1 + Afi) and N b s is the observed 
number density of galaxies. Since we assume that the 
clustering and the lensing are both weak, the clus- 
tering term A[i and the lensing term (5s — 2)k are 
additive and we get an observed correlation function 

, , / N obs (9 )N obs (9' + 9) \ 

w{0) 55 \ ni 7 

= w aa (6) + (5s - 2) 2 w RK (6) + 
2(5 S -2)uj gK (6). (10) 
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The brackets denote a directional average. It can be 
seen that there are three terms to the observed corre- 
lation function. The first term is due to the true clus- 
tering of galaxies. The second term is due to the mass 
density fluctuations while the last is a cross term. For 
details of the calculation of the effects of gravitational 
lensing, see V96. 

In order to evaluate the importance of magnifica- 
tion bias, we compare the two terms in Eq. |^ on a 
scale of 5 arcmin. The rms density fluctuation on that 
scale is approximately 0.035 for the sample R < 25.5 
(BSM). If the redshift of that sample is z ~ 1 then the 
rms k smoothed on the same scale is a K rs 0.015r2o cr 8- 
If s = 0.3 then the lensing term in Eq. ^ will be ap- 
proximate (5s — 2)k w 0.008 for as^o = 1- This is not 
negligible compared to the intrinsic clustering term 
unless f2o<7g <C 1. If our galaxy sample has a number- 
count slope of s = 0.2, we greatly increase the effect 
of lensing. In that case, the lensing term (5s — 2)« 
will double to w 0.016. Therefore, w(0) can be signif- 
icantly influenced by magnification bias as shown in 
fig. 1 of V96. In that paper the effect was evaluated 
in terms of median redshift of the sample. 

For a given set of cosmological parameters, such 
as Qq, A, us, and intrinsic galaxy-galaxy £, gg (r, z) and 
mass-mass £ mm (r, z) correlation functions we can cal- 
culate oj gg (9) and (5s — 2) 2 u KKl {0) for an observed 
redshift distribution N(z) of galaxies. In order to 
calculate 2 (5s — 2)uj gK {9), it is necessary to compute 
the galaxy-mass £ 9m (r, z) correlation function, or in 
other words, the relation between the galaxy and mass 
distributions. If there is no correlation, then uj gK is 
by definition zero. The next simplest assumption is 
the linear bias model, whereby the galaxy overden- 
sity equals the mass overdensity apart from a factor 
b. This is the model we use for theoretical predictions, 

CgmiXi ^0 — ^ £mm(r, j ^ggi^i Z^) — b £771771(^*5 ^0? 

(11) 

where b is the biasing factor. In the linear bias model 
we then get (V96) 

N(z) + 3 % [5s - 2] <j>{z) y(z) [1 + z] ff" 1 ^)] 
o 

x x(z) 1 "' (l + z) 7 ^ e ). (12) 

This is a generalization of equation^. Here x(z) and 
y(z) are the comoving radial and angular diameter 



distances and <f>{z) is the lensing selection function 
which is the integral over all sources more distant than 
z of the ratio of the comoving angular lens-source dis- 
tance Uls{z', z) and observer-source distance yos( z ) 

poo I I \ 

4>{z')= / dzN(z) yLS[Z ' Z) . (13) 

Jz> yos[z) 

As stated in V96, the cosmological constant A 
has little influence on weak gravitational lensing and 
therefore magnification bias. In Eq. [l2| the depen- 
dence of w{6) on the rms mass density fluctuations 
parameterized by a$ has been expressed in terms of 
the galaxy correlation length ro and the bias factor 
b for a particular galaxy population, using the linear 
bias model. 

The contribution to w(0) from gravitational lensing 
depends on the product O x tr 8 . For the cross term, 
the dependence is linear and for the pure lensing term 
the dependence is quadratic. In this paper we have 
performed all calculations for a model in which Qq = 
1 ; A = and for various values of as- We note that 
for a fixed value of ro , a scaling in as is equivalent to 
a scaling in b. 

In Figure 3 we show the theoretical predictions for 
the correlation amplitude including the effects of mag- 
nification bias. The slope of the number counts in the 
R band is measured to be s = 0.31 ± 0.03, consistent 
with results from brighter samples (Smail et al. 1995). 
The observed data points are shown in the two panels 
of Figure 3 together with the theoretical predictions 
both with and without magnification bias. In the left 
panel we assume that e = and show the predictions 
for the amplitude for er 8 = 0.5, i.e. a highly biased 
model (long-dashed curve) and for cr 8 = 1 (short- 
dashed curve), i.e. a unbiased model. In the right 
panel it is assumed that e = 0.8. 

These theoretical curves can be compared with the 
curves in V96, Figure 1. In that paper, the effects 
of magnification bias are shown in terms of median 
redshift, while here we show the effects in terms of 
limiting magnitude which is an observable. 

The basic effect of magnification bias is to decrease 
the observed amplitude of clustering. It is larger for 
small values of ro, and larger values of e. The effects 
of magnification bias is also an increasing function of 
limiting magnitude. These effects can be simply un- 
derstood in terms of a competition between the effects 
of true clustering and the clustering induced by the 
magnification bias. The intrinsic clustering contribu- 
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tion decreases with decreasing value of r and increas- 
ing values of e and limiting magnitude. The clustering 
contribution from magnification bias is an increasing 
function of <r$ and limiting magnitude. From the fig- 
ure, we find that the decrease in amplitude is in the 
range of 10% to 25% and is not very sensitive to the 
limiting magnitude. The amplitude is already high at 
a magnitude limit R ~ 25 which corresponds to the 
spectroscopic limit of large telescopes. As discussed 
below, this fact provides us with a possible way of 
directly measuring the magnification bias and there- 
fore f^oos in deep redshift surveys. Even though the 
effects of magnification bias are not negligible, they 
are not large enough to change the conclusions in sec- 
tion 4. 

In our analysis we have considered only the effect 
of magnification bias caused by large scale structure. 
Weak lensing by individual galaxy halos, i.e. galaxy- 
galaxy lensing (eg. Brainerd, Blandford, Smail 1996) 
can change the correlation function on small angular 
scales through the magnification bias. Given the un- 
certainty in the size and masses of galaxy halos, the 
magnitude of this effect is uncertain. However, judg- 
ing by the results of Brainerd, Blandford, & Smail 
(1996), which reported a weak detection of the effect, 
it is unlikely that this effect will make a significant 
contribution to w(9) . 

The reason the effects of magnification bias are not 
more noticeable is that there is a large fraction of 
low redshift galaxies in the sample for which magni- 
fication bias is not important. Furthermore, the full 
sample considered so far has a number-count slope 
s » 0.3, close to the neutral slope s = 0.4. This 
motivates us to look for a subsample for which the 
number-count slope is significantly lower than 0.3. 

6. Magnification Bias for a Color-selected 
Sample 

It is clear from Eq. [To] that the effects of magnifica- 
tion bias are larger if the number counts slope for the 
sample considered, s, is significantly different from 
0.4 . In particular V96 has shown that for a sample 
with s = 0.2, magnification bias can produce a dis- 
tinctive upturn in the amplitude of w(8) for samples 
with median rcdshifts z > 1. If the median redshift 
increases with limiting magnitude of the sample, then 
we can in principle expect to sec this distinctive up- 
turn at faint magnitudes. 

It has been noted by Broadhurst (1996) that the 



reddest galaxies, as defined by V — I colors, have 
a small number-count slope of s < 0.2. Moreover 
Broadhurst et al. (1996) have shown that the num- 
ber counts slope is a decreasing function of V — I 
color, i.e. red-selected samples have a shallower slope. 
Bearing this in mind, we define color selected subsam- 
ples such that the slope of the number counts is 0.2 . 
This was done as follows. For a given magnitude limit 
in R we arranged the galaxies in order of increasing 
R — I color index. We removed blue galaxies until the 
slope of the number counts of the remaining sample 
is s = 0.2 and calculated w(0) in the same way as for 
the full sample. 

The results are shown in Figure 4 in the same form 
as in Figure 1 , and are listed in columns 6 to 8 of Ta- 
ble 1. As for the full sample, we detect a correlation 
signal down to R ~ 29. Figure 5 is equivalent to fig- 
ure 3 for the full sample except that there are no BSM 
data points as they have not color-selected samples. 
The dotted curves are the same for the full sample 
shown in Figure 3, as they do not include the mag- 
nification bias. We have assumed the same redshift 
distribution for the color selected sample for illustra- 
tive purposes only. The two samples need not have 
the same redshift distribution. The dashed curves in- 
clude the effects of magnification bias. 

The theoretical predictions for the correlation am- 
plitude show the expected larger influence of the mag- 
nification bias. The major differences in including the 
magnification bias for the red-selected sample are: 1) 
At brighter magnitudes, R 25, the effect is signifi- 
cantly larger than for the full sample, but has a similar 
behavior. 2) At faint magnitudes the effects are much 
more complicated. For as = 0.5 models, the observed 
effect of magnification bias for the color selected sam- 
ple is larger than for the full sample, but again not 
enough to significantly affect the conclusions about 
To for a given e. For as = 1 models, the behavior 
is quite different with the curves flattening out and 
in some cases showing an upturn at faint magnitudes. 
Unfortunately, the upturn is not a generic feature but 
occurs only for models with weak intrinsic clustering. 
Furthermore, an upturn in the correlation amplitude 
at faint magnitudes could have other explanations, 
such as a very faint local population of galaxies. 

It is important to note that the effects of magnifi- 
cation bias can be larger at brighter magnitudes than 
at fainter magnitudes. This is important as it allows 
us to use spectroscopic information to measure the 
magnification bias. There are several ways of implc- 
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meriting this. If we have some redshift information, ei- 
ther photometric or spectroscopic, physical pairs can 
be removed from the sample and the observed w{6) 
would be a direct measure of the magnification bias. 
An alternative way is having spectroscopic redshifts 
for a random subsample so that we can derive N(z) 
and £ (r, z) directly. From this information we can 
compute the true w(G) from Limber's equation and 
compare this result with the measured w(9). The 
difference would then be the contribution from mag- 
nification bias from which we can measure fl^as- 

7. Conclusions 

Analysis of the multicolor HDF data using the two- 
point angular correlation function leads us to the fol- 
lowing conclusions: 

1- We detect a clustering signal down to a magnitude 
limit of R = 29. 

2- We find values of the amplitude of the 
w{6) consistent with those obtained from ground- 
based observations at R ~ 26. The amplitude contin- 
ues to decrease down to the faintest magnitude limit 
considered with a slope consistent with 7 ~ 1.8. 

3- Our results show that the measured amplitudes of 
w{9) are consistent with some of the theoretical mod- 
els. The best fit model is one with linear evolution 
of the correlation function i.e. e = 0.8 and a present 
day correlation length r ~ 4 h^ 1 Mpc similar to that 
observed for IRAS galaxies. This would be consistent 
with a large population of normal galaxies already 
being in place at z ^ 1.7 . 

4- We also show that even though magnification bias 
is expected to have an impact on the angular corre- 
lation function at faint magnitudes, it is not strong 
enough to seriously affect the the conclusions about 
r and e, even for a red selected sample, for which 
magnification bias should be more important. 

5- The effects of magnification bias can be important 
for samples with a limiting magnitude of R ^ 25, 
which corresponds to the spectroscopic limit of large 
telescopes. This makes it possible to measure the ef- 
fects of magnification bias and therefore get an alter- 
native way of measuring f^ocs- Given the estimates 
for the magnification bias, this requires a galaxy sam- 
ple large enough to allow measurements of w(9) to 
better than 10% accuracy. 

Deep imaging and redshift surveys from ground 
based telescopes will make it possible to draw more 
definite conclusions about the effects of magnification 



bias and its potential use for the measurement of 0. 
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Table 1: Observed correlation functions 
Mag Limit Rcdshift # of Galaxies Amplitude Error of Galaxies Amplitude Error 

Estimate Full Sample RedSample 



20.0 0.25 

21.0 0.35 

22.0 0.48 

23.0 0.64 

24.0 0.87 

25.0 1.11 



26.0 


1.35 


263 


0.19 


0.18 


249 


0.35 


0.19 


26.5 




384 


0.31 


0.14 


362 


0.30 


0.15 


27.0 


1.54 


536 


0.06 


0.10 


510 


0.10 


0.10 


27.5 




692 


0.15 


0.07 


655 


0.08 


0.08 


28.0 


1.71 


900 


0.13 


0.06 


839 


0.13 


0.06 


28.5 




1230 


0.071 


0.043 


1082 


0.073 


0.049 


29.0 


1.87 


1559 


0.063 


0.033 


1226 


0.057 


0.043 


29.5 




1732 


0.059 


0.030 


1256 


0.043 


0.041 



Figure Captions 

Figure 1. Measured w(8) for all galaxies brighter 
than the limiting magnitude R as indicated in each 
panel. Error bars represent the la Poisson errors. 
The solid lines are the best fits for fixed power law 
index 7 = 1.8. The integral constraint (see text) has 
been included in the fits. 

Figure 2. The observed amplitude of w(6), mea- 
sured at 1 arcsec separation as a function of limiting 
R magnitude, for 18 < R < 29. Filled circles are 
from HDF data and open circles are adapted from 
Figure 2 of BSM, including 1 a error bars are shown. 
The dashed line is a power-law fit to the data from 
BSM. The curves are theoretical predictions assuming 
N(z) as described in the text for different values of 
the present correlation lengths ro as indicated, with 
the amplitude increasing with the adopted ro. The 
dotted curves are for e = 0, while the solid curves are 
for e = 0.8 . 

Figure 3. Observed amplitude of w(9), measured 
at 1 arcsec separation as a function of limiting R mag- 
nitude, for 24 < R < 29. Points with filled circles 
are from HDF data and open circles are from BSM 
data. The curves are theoretical predictions assuming 
N(z) as described in the text for different values of 
the present correlation lengths ro as indicated. In the 
left panel e = is assumed while in the right panel 
e = 0.8 is assumed. The dotted curves are predictions 



without the effects of magnification bias. The dashed 
curves are equivalent curves including the effects of 
magnification bias. The long dashed curves are for 
as = 0.5 and the short dashed curves are for as = 1. 
The observed number counts slope s = 0.3 is used in 
the predictions. 

Figure 4. Measured w(9) for the red selected 
galaxies brighter than the limiting magnitude R as 
indicated in each panel. Error bars represent the la 
Poisson errors. The solid lines are the best fits for 
fixed power law index 7 = 1.8. The integral con- 
straint (see text) has been subtracted from the fits. 

Figure 5. Observed amplitude of w(9), measured 
at 1 arcsec separation as a function of limiting R mag- 
nitude, for 24 < R < 29 for the red-selected sam- 
ple. The data points are from the HDF data sample 
only. The curves are theoretical predictions assuming 
N(z) as described in the text for different values of 
the present correlation lengths ro as indicated. In the 
left panel e = is assumed while in the right panel 
e = 0.8 is assumed. The dotted curves are predictions 
without the effects of magnification bias. The dashed 
curves are equivalent curves including the effects of 
magnification bias. The long dashed curves are for 
a 8 = 0.5 and the short dashed curves are for as = 1. 
The observed number counts slope s = 0.2 is used in 
the predictions. 



10 



